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ABSTRACT 


Accurate  physical  scale  models  of  complex  systems  are  found  to 
be  cf  great  use  infields  as  diverse  as  aerodynamics,  fluid  mechanics, 
radar,  or  nuclear  physics.  It  now  appears  that  the  concept  of  scale 
modeling  may  offer  distinct  advantages  to  those  interested  in  the  study 
of  thermodynamic  processes  that  occur  in  large  physical  structures.  In 
this  investigation,  it  is  shown  from  a  study  of  the  heat  equation  and  its 
boundary  conditions  that  physical  scale  modeling  can  be  used  to 
simulate  realistic  systems  operating  in  realistic  outdoor  environments. 
The  thermal  properties  of  the  construction  materials  used  in  the  system 
under  study  are  allowed  to  vary  with  position,  thus  allowing  the 
structure  to  be  divided  into  a  number  of  compartments.  Physical  effects 
which  involve  heat  exchange  between  the  structure  and  the  ocean,  or 
with  the  atmosphere,  are  investigated.  Both  time-dependent  and  time- 
independent  cases  are  examined.  The  scaling  laws  of  thermodynamic 
similitude  which  govern  the  exchange  of  infrared  radiation  between  any 
structure  and  its  environment  are  derived.  These  scaling  laws  are  then 
used  to  obtain  a  set  of  dimensionless  variables  which  transform  the  heat 
equation  and  its  boundary  conditions  into  manifestly  scale  invariant 
form.  In  this  form  the  heat  equation  applies  to  a  scale  model  of 
arbitrary  size.  Consideration  of  the  "temperature"  Junction  solution  to 
the  scale  invariant  heat  equation  is  applied  to  the  design  of  an  actual 
physical  scale  model  and  to  the  simulation  of  a  background  environment 
in  which  the  thermal  properties  of  die  model  are  to  be  measured. 
Methods  for  the  calibration  and  validation  of  laboratory  measurements 
using  actual  field  measurements  are  described.  The  scaling  laws  that 
are  derived  apply  to  realistic  physical  structures  interacting  with  an 
outdoor  environment.  In  particular,  they  may  be  used  to  accurately 
describe  such  structures  as  aircraft  inflight,  ships  at  sea,  tanks,  or 
even  buildings. 


ADMINISTRATIVE  INFORMATION 

This  work  is  in  partial  fulfillment  of  Milestone  2,  Task  3,  of  the  Topside 
Signature  Reduction  Project  (RH21C17)  of  die  Surface  Ship  Technology  Program 
(SC1A/PE0602121N).  The  work  described  herein  was  sponsored  by  the  Office  of  the 
Chief  of  Naval  Research  (OCNR  33)  and  was  performed  by  the  Carderock  Division, 
Naval  Surface  Warfare  Center  (CDNSWC,  Code  7230). 
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INTRODUCTION 


The  purpose  of  this  investigation  is  to  obtain  a  better  understanding  of  how 
thermal  radiation  is  exchanged  between  a  physical  structure  and  its  environment.  At 
presort,  such  studies  are  performed  using  experimental  data  and  calculated  results 
based  upon  computer  codes.  It  is  not  generally  recognized  that  valuable  information 
can  be  obtained  from  laboratory  measurements  made  on  scale  models.  Much  can  be 
learned  from  the  uses  to  which  scale  models  have  long  been  put  in  other  fields  of 
study;  among  others,  in  the  areas  of  radar,  acoustics,  fluid  mechanics,  and 
aerodynamics.  In  the  field  of  thermal  radiation  as  well,  it  is  thought  that  scale  models 
can  be  a  valuable  investigation  resource.  The  use  of  scale  models  has  the  potential  for 
cost  effectively  enhancing  the  value  of  field  measurements  and  computer  code 
calculations.  For  example,  scale-model  studies  may  be  used  to  verify  those  carried  out 
in  the  field  and  provide  important  consistency  checks  with  which  computer  codes 
must  comply.  This  should  lead  to  improvements  in  present-day  computer  codes  based 
upon  a  more  complete  representation  of  the  physical  laws  that  give  rise  to  the 
"temperature"  function,  which  ultimately  is  responsible  for  the  way  a  physical  system 
radiates  in  the  infrared. 

The  feasibility  of  scale  modeling  is  examined  theoretically  by  analyzing  the 
heat  equation.  This  equation,  together  with  the  appropriate  boundary  conditions, 
offers  a  physically  accurate  description  of  all  the  convective  and  radiative  processes 
that  take  place  whatever  a  physical  structure  exchanges  heat  with  its  surroundings. 

For  instance,  one  such  physical  structure  may  be  a  ship  sailing  through  the  ocean,  an 
airplane  flying  through  the  air,  or  a  tank  moving  through  a  smoke-charged  wind.  The 
equations  are  written  to  be  applicable  to  thermodynamic  systems  that  include 
compartments  made  from  materials  that  may  vary  so  that  the  thermal  parameters  that 
characterize  these  materials  can  also  vary.  The  conclusions  from  this  investigation  of 
the  heat  equation  cover  the  time-dependent  as  well  as  the  time-independent  cases. 

The  results  of  analyzing  the  heat  equation  are  the  thermodynamic  scaling  laws 
appropriate  for  thermodynamic  scaling.  In  the  case  of  a  ship,  for  example,  the  scaling 
laws  enable  experimental  results  from  scale  models  to  be  compared  with  those  from 
full-size  ships  at  sea.  The  scaling  of  temperature,  position,  and  time  is  related  to 
scaling  of  heat  sources,  boundary  conditions  (including  initial  conditions),  and 
material  thermodynamic  parameters.  Thus,  true  scaling  of  a  large  ship  involves  much 
more  tten  simply  bu:lding  a  small  geometrical  replicate.  Thermodynamic  similitude 
requires  that  response  times,  heat  loads,  thermal  properties  of  ship  materials  and  of 
fluids  in  the  ship’s  surroundings  be  taken  into  account  as  well. 

Furthermore,  the  scaling  laws  suggest  a  set  of  dimensionless  variables  that 
transform  the  heat  equation  into  a  manifestly  scale  invariant  form,  which  describes  all 
geometrically  and  thermodynamically  similar  ships  no  matter  what  their  sizes.  Such  a 
description  is  hence  size-independent.  The  set  of  dimensionless  variables  define  the 
minimum  number  of  measurements  that  fix  the  "temperature  function"  and,  therefore, 
define  the  most  cost-effective  procedure  for  the  use  of  scale  models.  The  set  of 
dimensionless  variables,  once  obtained,  are  in  themselves  implementation  criteria 
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serving  to  support  a  wide  variety  of  applications  which  may  be  envisioned,  such  as 
the  setting  of  thermal  radiation  goals  for  surface  ships,  the  study  of  anti-ship  seekers, 
engineering  design  choices  related  to  ship  size,  and  physical  concepts  that  should  lead 
to  the  next  generation  of  thermal  computer  programs. 


THERMAL  (INFRARED)  SCALING  LAWS: 

THE  TIME-DEPENDENT  CASE 

Scaling  laws  appropriate  to  a  structure  based  on  compartments  made  from 
materials  with  thermal  properties  that  are  allowed  to  vary  from  place  to  place  are 
derived.  The  compartments  are  filled  with  fluids  (air,  water,  oil,  etc.).  Within  the 
compartments  both  convective  and  radiative  heat  transfer  processes  occur.  All 
conductive  processes  that  occur  throughout  the  structure  are  described  by  the  heat 
equation.  Equation  (1).  Of  all  possible  solutions  to  the  heat  equation,  those  applicable 
to  a  particular  structure,  in  a  particular  environment,  are  selected  by  the  prevailing 
boundary  conditions.  The  heat  transfer  processes  between  compartments  and  their 
enclosing  conductive  walls  are  accounted  for  by  the  boundary  conditions  that  are 
imposed.  11131  is,  the  fluid-filled  space  of  compartments  are,  from  the  heat  equation 
describing  the  heat  conduction  processes  in  such  a  system,  part  of  the  internal 
thermodynamic  surroundings.  The  fluids  on  the  outside  (the  air  of  the  atmosphere 
and  the  water  of  the  ocean)  compose  the  external  thermodynamic  surroundings.  The 
total  thermodynamic  surroundings  (the  internal  and  external  portions  together)  are 
separated  from  the  conductive  thermodynamic  system  by  a  surface  termed  the 
thermodynamic  boundary.  The  heat  exchange  processes  that  occur  connecting  the 
thermodynamic  system  to  any  portion  of  its  surroundings  are  described  by  means  of 
appropriate  boundary  conditions.  The  heat  equation  contains,  as  parameters,  the 
thermodynamic  properties  of  materials  from  which  the  structure  is  made.  To  be  as 
realistic  as  possible,  it  is  necessary  to  consider  these  thermodynamic  properties  to  be 
completely  general  functions  of  position. 

The  mathematical  structure  containing  the  thermodynamic  description  above  is 
given  by  a  form  of  the  heat  equation  and  boundary  conditions  that  can  be  displayed 
as: 


V  •[*(*) V7fc,f) ]  ♦  HQc,f)  =  p(x) c(x)  ^-(x,/)  (1) 

k(x)  =  /j(x)  AT\x,t)  for  all  x  C  b  (2) 

for  all  x  0) 
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The  first  equation  is  the  heat  equation,  the  second  is  a  boundary  condition,  and  the 
third  is  an  initial  condition.  The  various  symbols  are  defined  as: 


boundary 

three-dimensional  gradient  operator 

operator  yielding  gradient  magnitude  in  direction  normal  to  boundary 
thermal  conductivity  (power/distance/temperature) 
temperature  at  position  x 

coordinates  of  a  point  in  three-dimensional  position  space 
coordinates  of  a  point  in  time 

rate  jf  heat  production  per  unit  volume  at  position,  x,  and  time,  t, 
(power/volume) 

mass  density  at  position  x  (mass/volume) 

specific  heat  at  position  x  (heat/temperature/mass) 

heat  transfer  coefficient  at  position,  x,  belonging  to  boundary,  b 
(power/area/temperature  differential) 

temperature  difference  prevailing  between  ship  and  its  surroundings  at 
position,  x,  and  time,  t 

initial  temperature  distribution  that  holds  at  every  position,  x, 
throughout  the  ship 


The  solution  of  Equations  (1)  to  (3)  yields  a  temperature  function 

T  *  r(r,  t,  H,  k,  p,  c,  ht  A T)  (4) 

which  evolves  from  its  initial  value  T0.  It  is  a  function  of  position  and  time  which 
depends  upon  the  heat  sources,  the  material  thermodynamic  properties,  the  heat 
transfer  coefficient,  and  the  temperature  differences  between  system  and  surroundings. 

The  material  thermodynamic  properties,  k,  p,  and  c,  are  completely  general 
functions  of  position,  as  is  the  heat  transfer  coefficient,  h.  Note  that  the  boundary 
condition,  Equation  (2),  is  Newton’s  law  of  heat  transfer  and  properly  describes 
natural  convection  and  forced  convection  to  fluid  surroundings  (air,  water,  oil,  etc.) 
as  well  as  radiation  in  accordance  with  Planck’s  frequency  distribution  or  Stefan’s 
law.  In  brief,  the  complex  processes  of  heat  transfer  between  a  ship  with  realistic 
architecture  and  materials  and  its  fluid  surroundings,  both  internal  and  external,  are 
accurately  described  by  the  heat  equation  and  its  companion  boundary  conditions. 
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The  scaling  laws  are  obtained  by  applying  the  heat  equation  and  boundary 
conditions  to  both  a  ship,  S,  and  its  physically  scaled  model,  M,  and  by  stipulating 
that  the  equations  hold  in  both  instances.  Geometric  similarity  of  model  and  ship  is 
obtained  by  enforcing  the  relation  between  analogous  points  x*1  and  xs  of  model  and 
ship 


(5) 


where  Kx  is  a  scaling  constant,  namely  that  of  position.  Complete  infrared  similarity 
of  model  and  ship  is  obtained  by  also  requiring  analogous  relations  that  involve  all 
pertinent  quantities  from  the  heat  equation  and  boundary  conditions.  Thus,  Equations 
(6)  through  (13)  define  the  scaling  constants  of  time,  temperature,  conductivity,  mass 
density,  specific  heat,  heat  transfer  coefficient,  temperature  difference,  and  initial 
temperature,  respectively. 


tM 

Kt  =  *—  (6) 

ts 

«l»Af 

K  T  *  (7) 

j '5 


(8) 


(9) 


(10) 


(ID 
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kat  = 


(12) 


atm 

ATS 


Equations  (1)  through  (3)  applied  to  a  ship  yield: 


Vs  (ksVsTs)+Hs  =  pscs— 

dts 


(13) 


(14) 


ksV^Ts  =  hsATs  for  all  xs  C  bs  (boundary  of  full-scale  ship)  (15) 


=  T0*  for  all  xs  (16) 

where  all  symbols  retain  their  previous  meanings.  Superscript  "S"  indicates  quantities 
pertaining  to  a  ship,  and  for  brevity  of  notation  space  and  time  dependence  is  not 
shown  explicitly  although  it  is  implicitly  understood.  Similarly,  for  an  analogous 
model, 


V*(kMVMTM)+HM  9  pMcMdT * 

dtM 


(17) 


kMV^ Tm  =  huATu  for  all  xM  C  bM  (boundary  of  scale  model)  (18) 


Tm =  T0m  for  all  xM  (19) 

By  invoking  the  scaling  constants,  Equations  (6)  through  (13),  the  ship 
equations  can  be  rewritten  as: 
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(20) 


^£rVM.(k»VuTM),KHHu  =  K“KfT 
Kx2  Kt 


K*KT],mvMtm  =  Kf^AJh^ATM  for  all  xMCbM  (21) 

Kr 


Kt  TM\tHxo  =  KtT0m  for  all  xm  (22) 

But  by  insisting  that  the  ship  equations,  (20)  to  (22),  and  the  model  equations,  (17)  to 
(19),  must  both  hold,  requires  that: 


K  k  — r  =  1  (23) 

—  K  2 
pc  "-x 


Wt  .  ! 

Vt 


(24) 


KhKx^  =  l  (25) 

*  KT 


(26) 


The  abbreviating  definitions  are  obvious  in  context. 

Equations  (23)  through  (26)  are  the  thermal  scaling  laws  for  the  time- 
dependent  case.  They  connect  the  temperatures  to  scaling  of  time  and  distance, 
conductivity,  density,  specific  heat,  heat  sources,  heat  transfer  coefficients, 
temperature  differences  between  boundaries  and  surroundings,  and  initial  temperature 
distributions.  These  laws  apply  to  realistic  naval  ship  architectures  and  materials  in 
die  context  of  heat  transfer  mechanisms  including  conduction,  convection,  and 
radiation.  They  also  are  the  basis  upon  which  an  IR  physical  scale  modeling 
measurements  program  can  be  designed. 
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DIMENSIONLESS  VARIABLES  ASSOCIATED  WITH  THERMAL 
(INFRARED)  PHYSICAL  SCALE  MODELING 

Essentially,  physical  scale  models  can  be  built  to  thermally  behave  like  full- 
scale  ships  because  the  thermodynamic  equations  can  be  cast  entirely  in  terms  of 
dimensionless  variables.  By  definition,  these  are  variables  with  magnitudes 
completely  independent  of  size,  and  they  are  exactly  the  same  applied  to  a  ship  and  to 
physical  scaled  models.  An  appropriate  set  of  dimensionless  variables  and  parameters 
is  immediately  suggested  by  the  form  of  the  scaling  laws  previously  derived.  From 
such  dimensionless  quantities,  the  entire  set  of  thermodynamic  equations  can  be  cast 
into  dimensionless  form.  Thus,  a  ship  and  its  scaled  model  behave  according  to 
equations  manifestly  independent  of  size.  This  leads  to  the  logical  design  of  a 
physical  scale  modeling  measurements  laboratory  because  the  physical  variables  that 
such  a  laboratory  must  control  and  measure  include  only  the  dimensionless  quantities 
themselves.  True  physical  scale  modeling  is  therefore  tantamount  to  a  laboratory 
arrangement  that  fixes  all  dimensionless  quantities  to  be  the  same  values  on  a  model 
and  the  ship  it  represents  by  virtue  of  displaying  precisely  the  same  IR  dynamics  as 
the  ship  itself. 

To  obtain  appropriate  dimensionless  quantities,  the  heat  equation  and  its 
boundary  conditions  are  cast  in  dimensionless  form.  The  dimensionless  quantities 
required  by  the  scaling  of  the  IR  problem  are: 


x  = 


x 

L 


(27) 


1  =  — — —  (Fourier  number)  (28) 

pc  l1 


T=L  (29) 

*o 


=  “!  (30) 

kTo 


h  =  —  (Nusselt  number)  (31) 

where  the  bar  symbol  indicates  a  dimensionless  quantity,  and  L  represents  a  "typical" 
length;  all  other  symbols  are  as  previously  defined.  The  five  dimensionless  variables 
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in  Equations  (27)  to  (31)  are  suggested  immediately  by  the  five  scaling  laws  derived, 
viz.,  the  condition  of  geometric  similarity.  Equation  (5),  and  the  four  conditions  of 
thermodynamic  similarity  which  follow  from  the  heat  equation  and  its  boundary 
conditions,  Equations  (23)  to  (26). 

Thus,  x,  which  is  consistent  with  Equation  (S),  is  simply  a  reformulation  of 
the  condition  of  geometric  similarity,  and  is  obviously  dimensionless.  In  T,  a 
definition  is  immediately  suggested  by  the  cluster  of  subscripts  in  Equation  (23).  That 
t  is,  in  fact,  dimensionless  follows  from  a  detailed  inspection  of  units  carried  by  the 
quantities  contained  in  its  definition.  T  is  suggested  by  the  cluster  of  subscripts  in 
Equation  (26),  and  also  is  obviously  dimensionless.  H  is  suggested  by  the  subscript- 
cluster  in  Equation  (24),  writing  t  and  T  in  terms  of  their  dimensionless  analogues. 
Thus, 


l  =  tpcL2  1 
pc  T  pc  k 


(32) 


which  simplifies  to 


l  _  HL?  t_ 
pc  T~  kTa  f 


(33) 


In  Equation  (33)  the  coefficient  of  dimensionless  t  /T  defines  H,  whose  own  form  is 
dimensionless  as  follows  from  a  detailed  inspection  of  units  carried  by  the  quantities 
entering  its  definition.  Similarly,  IT  is  suggested  by  the  cluster  of  subscripts  contained 
in  Equation  (25),  from  which: 


(34) 

k  T  k  t  t 
‘o1 

and  on  simplifying, 

h  T  AT  _  hLAT 

k  T~  k  f  ™ 

Here  IT  as  appears  as  the  coefficient  of  dimensionless  AT/T.  Inspection  indicates  If 
also  is  dimensionless,  as  required.  Thus,  the  form  of  the  IR  scaling  laws  suggests  the 
form  of  the  dimensionless  quantities,  which  are  appropriate  to  the  scaling  of  IR. 
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All  of  the  dimensionless  quantities  in  Equations  (27)  through  (31),  have  well- 
defined  physical  interpretations.  The  Fourier  number  measures  the  ratio  of  heat 
conducted  through  a  volume  element  to  that  "stored,”  causing  an  increase  in 
temperature.  The  Nusselt  number  ¥  measures  the  ratio  of  heat  "convected”  from  a 
surface  to  that  conducted  through  the  surface.  Less  recognizable,  H  for  an  elemental 
volume  is  the  ratio  of  heat  from  a  source  present  to  that  "stored,”  causing  an  increase 
in  temperature.  More  obviously,  IT  is  a  measure  of  position  in  units  of  a 
representative  length,  and  T  of  temperature  in  units  of  the  initial  temperature. 

Using  the  dimensionless  quantities  suggested  by  the  scaling  laws  naturally 
transforms  the  thermodynamic  description  of  a  ship  interacting  with  its  surroundings 
into  one  that  is  dimensionless,  i.e.,  completely  independent  of  size  or  scale.  Using 
the  relations 


and 


d_  m  _k_  d_ 
dt  pcL 2  ft 


(36) 


(37) 


which  follow  from  Equations  (27)  and  (28)  and  replacing  T,  H,  and  h  with  their 
equivalents  in  terms  of  T,  H,  and  ¥  given  by  Equations  (29)  to  (31),  Equations  (1)  to 
(3)  become: 


V2  T  +  —  V  T  +  H  =  LI  (38) 

*  dl 


VnT  =  hAT  for  all  x  C  b  (39) 

and 

T\~tmo  =  1  for  all  x  (40) 

Equations  (38)  to  (40)  are  the  full  information  equivalents  of  the  heat  equation  and 
boundary  conditions,  but  are  now  in  manifestly  scale  invariant  form.  Every  term  is 
dimensionless,  thus  corroborating  the  choice  of  dimensionless  variables  and 
parameters,  Equations  (27)  to  (31),  suggested  by  the  form  of  the  IR  scaling  laws. 
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The  scale  invariant  heat  equation  and  boundary  conditions  imply  the  existence  of  a 
"temperature"  function, 

T  =  T  (x,  7,  H,  h,  AT)  (41) 

which  evolves  from  its  initial  value  of  unity.  It  is  a  function  of  "position"  and  "time" 
which  depends  parametrically  upon  the  "heat  sources,"  the  "heat  transfer  coefficient," 
and  the  "temperature  differences"  between  systems  and  surroundings.  This 
expression  for  the  "temperature"  is  the  analogue  of  the  temperature  solution  given  by 
Equation  (4),  but  here  given  in  manifestly  scale  invariant  form.  Moreover,  x  displays 
an  efficient  reduction  in  the  number  parameters  upon  which  it  depends,  when 
compared  to  the  case  of  T. 

The  temperature  as  given  by  Equation  (4)  depends  upon  eight  parameters,  viz., 
x,  t,  H,  k,  p,  c,  h,  and  AT;  whereas,  the  scale  in  variant_"  temperature"  is  seen  to 
depend  upon  only  five  scale  invariant  parameters,  viz.,  x,  t ,  H,  TT,  and  AT.  The 
matter  of  reducing  the  number  of  parameters  from  eight  to  five  is  of  importance  to 
the  feasibility  of  any  measurements  program  which  seeks  to  effectively  study, 
experimentally,  the  solutions  of  the  heat  equation. 

Reviewing  the  development  thus  far  reveals  that  an  IR  physi'  A  scale  modeling 
program  for  ships  is  a  feasible  undertaking.  A  ship’s  temperature  function,  together 
with  the  emissive/ reflective  properties  of  the  outer  sldn,  ultimately  determines, 
through  Planck  radiation,  the  IR  signal  emitted  by  a  ship.  Thus,  the  problem  of 
studying  the  IR  emission  of  a  ship  is  reducible  to  obtaining  the  temperature  function 
of  the  ship.  But  this  is  given  as  a  solution  of  the  heat  equation  for  the  ship  together 
with  its  companion  boundary  conditions.  In  order  to  reduce  the  problem  of  studying 
IR  ship  emissions  to  laboratory  manageability,  the  heat  equation  must  equally  apply  to 
a  ship  and  its  geometrically  similar  physical  scale  model.  The  conditions  that  a  point- 
by-point  comparison,  in  space  and  time,  between  a  ship  and  model  shall  hold  are 
embodied  in  the  IR  scaling  laws.  But  once  a  model  is  built  in  accordance  with  the 
scaling  laws  another  logical  question  arises;  viz.,  what  is  the  most  efficient  set  of 
measurements  for  the  study  of  the  IR  defining  temperature  function?  The  answer  is 
obtained  by  casting  the  heat  equation  and  its  boundary  conditions  into  dimensionless 
form.  The  solution  is  a  temperature  function  parametrically  dependent  upon  a 
minimal  number  of  dimensionless  quantities.  In  this  way,  the  laboratory  effort 
required  to  study  the  temperature  function  is  optimized.  Physical  scale  modeling  in 
the  IR  becomes,  therefore,  a  laboratory  arrangement  to  equate  the  appropriate 
dimensionless  quantities  of  a  scale  model  to  those  of  the  actual  ship  under  study. 

Thus,  physical  scale  modeling  is  seen  to  be  an  activity  which  is  entirely  feasible  in 
accordance  with  standard  thermodynamic  practice,  and  holds  promise  of  being  an 
efficient,  cost-effective  means  of  studying  the  IR  contrast. 
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The  following  discussion  examines  particular  cases  in  which  certain 
dimensionless  variables  are  associated  with  convection  and  radiation  exchange  at  a 
boundary. 

True  scale  modeling  implies  that  the  Nusselt  number,  F,  defined  in  Equation 
(31),  shall  be  fixed  at  the  same  value  in  both  model  and  ship.  The  magnitude  of  IT 
controls  the  heat  exchange  mechanisms  everywhere  on  the  boundary,  as  may  be  seen 
from  the  boundary  conditions,  Equation  (39).  The  temperature  function,  T,  is 
parametrically  dependent  on  F,  as  indicated  in  Equation  (41). 

The  two  principal  mechanisms  of  heat  exchange  at  a  ship  boundary  are  those 
of  convection  and  radiation.  Convection  occurs  between  the  ship  and  the  sea,  and 
between  the  ship  and  the  air/atmosphere.  Depending  upon  operating  conditions  ar 
environmental  conditions,  free  convection  or  forced  convection  may  apply,  with  ei. 
laminar  or  turbulent  flow  conditions.  Radiation  exchange  occurs  between  the  ship  and 
the  air/atmosphere  and  between  the  ship  and  the  sun.  All  of  the  heat  exchange 
mechanisms  can  be  occurring  simultaneously  at  a  given  position,  or  separately  at 
different  segments  of  a  ship.  Thus,  all  of  the  mechanisms  must  be  under  the 
experimental  control  of  the  scale  model  laboratory.  The  Nusselt  number  F  can  be 
written  as  a  superposition  of  convection  and  radiation  terms.  Hence, 

h  -hc  +  hr  (42) 

It  can  be  established  purely  by  dimensional  analysis  that  the  convective  contribution  to 
the  Nusselt  number  has  a  functional  dependence  given  by 


hc  =  hc  (Re,  Pr,  Gr)  (43) 

where  the  variables  are  respectively,  the  Reynolds  number,  the  Prandtl  number,  and 
the  Grashof  number.  Each  of  these  is  a  dimensionless  number  and,  therefore,  Fc 
embodies  convective  effects  in  a  scale  invariant  form.  The  dimensionless  quantities 
are  defined  as: 


Re  -  v  (Reynolds  number)  (44) 

n 

Pr  *  IS  (Prandtl  number)  (45) 
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(Grashof  number) 


(46) 


Gf  -  fig  A  T  L  y 

n2 


where: 

P 

fluid  mass  density 

L 

= 

representative  length 

V 

= 

fluid  ffee-stream  velocity 

7 

= 

fluid  velocity 

c 

= 

fluid  heat  capacity 

k 

* 

thermal  conductivity 

0 

= 

thermal  expansion  coefficient 

g 

= 

acceleration  constant  due  to  gravity 

AT 

= 

temperature  differential 

Each  of  the  dimensionless  numbers  defined  in  Equations  (43)  through  (46)  has 

a  well-defined  physical  significance.  The  convective  portion  of  the  Nusselt  number, 
TTC,  is  a  measure  of  the  ratio  of  the  flux  of  power  convected  to  that  which  would  be 
conducted  were  the  fluid  at  rest.  The  Reynolds  number.  Re,  which  fixes  the  flow 
pattern  around  a  solid  body  in  a  moving  fluid,  measures  the  ratio  of  the  flux  of  fluid 
momentum  to  the  flux  of  fluid  viscous  drag.  The  Prandtl  number,  Pr,  relates  the 
temperature  distribution  in  a  fluid  to  its  velocity  distribution,  and  is  a  measure  of  the 
ratio  of  the  kinematic  viscosity,  ij/p>  to  the  thermal  diffusivity,  k/pc.  The  Grashof 
number,  Gr,  fixes  the  relative  importance  of  buoyancy  and  viscous  forces,  and  is  a 
measure  of  the  ratio  of  fluid  buoyancy  force  per  area  to  viscous  drag  per  area.  In  the 
case  of  pure  free  convection,  the  Reynolds  number  is  inconsequential,  and  the 
convective  Nusselt  number,  Equation  (43),  simplifies  to 


TTc  =  7Tc(Pr,  Gr)  (47) 

Similarly,  for  pure  forced  convection,  the  Grashof  number  drops  out  of  consideration 
yielding 


Tc  =  JTC  (Re,  Pr)  (48) 

This  last  expression  holds  for  either  laminar  or  turbulent  flow,  although  the  specific 
functional  dependence  will  change  between  die  two  types  of  flow.  If  the  fluid  to 
which  this  last  pair  of  equations  is  applied  is  a  gas,  a  further  simplification  occurs. 
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Over  a  wide  range  of  temperature,  the  Prandtl  number  is  practically  the  same  for  any 
gas  and  therefore  may  be  neglected,  obtaining  from  Equations  (47)  and  (48): 


K  -  K  (Gr)  (49) 

and 

rc  =  JTC  (Re)  (50) 


for  gases  only. 

The  radiation  contribution  to  the  Nusselt  number  has  the  functional  dependence 


K  =  TTr  (St)  (51) 

where  the  variable  is  the  Stefan  number  that  is  dimensionless  and,  therefore,  TTr 
embodies  radiative  effects  in  a  scale  invariant  form.  The  Stefan  number  is  defined  as: 


(Stefan  number) 


(52) 


where: 

or  =  Stefan  constant  from  the  T4  radiation  law 

€  *  emissivity 

T  =  temperature 

L  =  representative  length 

k  =  thermal  conductivity 

The  physical  significance  of  the  Stefan  number  is  that  it  measures  the  ratio  of  heat 
radiated  to  that  conducted  per  unit  area. 

As  stated,  true  physical  scale  modeling  requires  equating  the  total  Nusselt 
number.  Equation  (31),  in  ship  and  model.  Convective  and  radiative  exchange  at  die 
boundaries  each  contribute  to  the  overall  Nusselt  number.  In  the  most  general  case, 
these  contributions  may  be  accounted  for  through  control  of  four  dimensionless 
quantities,  viz.,  die  Reynolds  number,  the  Prandtl  number,  the  Grashof  number,  and 
die  Stefan  number,  of  Equations  (43)  and  (51).  In  the  particular  limiting  cases 
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discussed,  simplification  occurs  requiring  control  over  fewer  dimensionless  quantities, 
as  in  Equations  (47)  to  (SO). 

In  addition  to  equating  the  overall  Nusselt  number,  IT,  applied  to  ship  and 
model,  the  relative  contributions  of  convective  and  radiative  effects  should  also  be 
represented  properly.  True  scale  modeling  therefore  entails  separately  equating  TTC  in 
ship  and  model  as  well  as  Fr. 


THERMAL  (INFRARED)  SCALING  LAWS: 

THE  TIME-INDEPENDENT  CASE 

To  derive  the  scaling  laws  appropriate  for  a  ship  in  a  thermodynamic  steady 
state,  the  formalism  and  definitions  used  in  the  time-dependent  case  are  retained.  The 
physical  nature  of  the  ship  and  its  energy  exchange  processes  with  its  surroundings 
are  as  previously  discussed,  except  that  changes  in  time  are  not  considered;  no 
transients  are  allowed.  In  particular,  the  ship  is  made  of  materials  whose  conductivity 
is  allowed  to  vary  from  place  to  place,  and  has  a  naval  architecture  providing  for  a 
variety  of  compartments  filled  with  fluids  (air,  water,  oil,  etc.). 

The  mathematical  structure  capturing  the  thermodynamic  description  in  this 
steady  state  case  is: 


V  [*(x)VrCt)]+H(x)  *  0  (53) 

kfryVJ  (x)  =  h$)AT  (x)  for  all  x  C  b  (54) 

The  solution  of  the  time-independent  heat  equation  and  its  companion  boundary 
condition,  Equations  (53)  and  (54),  respectively,  yields  a  temperature  function 

T  =  T(x,  H,  k,  h ,  AT)  (55) 


which  is  a  function  of  position  only,  and  depends  parametrically  upon  the  heat 
sources,  the  thermodynamic  conductivity,  the  heat  transfer  coefficient,  and  the 
temperature  differences  between  system  and  surroundings.  Note  that  k  and  h  are 
general  functions  of  position,  and  Newton’s  law  again  describes  heat  transfer 
mechanisms  at  the  boundaries,  e.g.,  convection  and  radiation. 

As  before,  the  time-independent  scaling  laws  are  extracted  from  the  heat 
equation  and  its  boundary  conditions  simply  by  applying  them  to  both  a  ship  and  its 
physically  scaled  model.  Hence,  for  a  ship: 
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v5(*vsr)  +  hs  =  o 


(56) 


ksvt J5  =  for  aU  xs  C bs 


(57) 


and,  for  its  analogous  model: 


VM-(kMVMTM)  +  HM  =  0 


(58) 


kMV*fTM  =  hMLTM  for  all  xM  C  bM 


(59) 


where  all  symbols  in  Equations  (56)  to  (84)  retain  their  previous  meanings. 

Invoking  the  scaling  relations,  Equations  (5)  through  (13),  the  ship  equations 
become: 


VM  (kMVuTM)  +  KuHm  =  0 
Kx2 


(60) 


and 


Wt  kMvMTM  m  Kh  hmatm  for  all  xmC  bM  (61) 


For  the  ship  equations,  (60)  and  (61)  and  the  model  Equations  (58)  and  (59)  to  both 
hold,  it  must  be  that 


**  KT 


=  1 


(62) 


Kh  Kx  K&r  =  1  (63) 

1  T 

These  are  important  scaling  laws  that  must  hold  in  the  steady  state  case.  They 
connect  the  temperatures  to  the  scaling  cl  distance,  conductivity,  heat  sources,  heat 
transfer  coefficients,  and  temperature  differences  between  boundaries  and 
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surroundings.  These  laws  are  the  basis  for  the  design  of  physical  scale  modeling 
measurements  in  the  steady  state,  and  they  supply  restrictions  to  which  computer 
models  must  conform. 

The  dimensionless  variables  appropriate  to  the  steady  state  case  are  especially 
easy  to  obtain.  They  are  already  defined  in  Equations  (27)  through  (31),  with  the 
exception  of  Equation  (28)  that  depends  explicitly  on  time.  These  equations  now 
become: 


J  ■  z 


T  •  -Z-  (65) 

V 


H  =  T3T-  (66) 

kT0> 


(67) 

Here  T0.  is  redefined  to  mean  an  average  background  temperature.  The  prime  in  the 
subscript  symbol  is  a  reminder  of  the  new  definition  imposed.  The  physical 
interpretations  previously  given  for  the  dimensionless  quantities,  Equations  (64)  to 
(67),  are  retained.  As  expected,  these  dimensionless  quantities  are  adequate  to 
transform  the  steady  state  thermodynamic  description  of  a  ship  exchanging  heat  with 
its  surroundings  to  one  which  is  stale  invariant. 

Using 


V  =  1  V  (68) 

L 4 

and  replacing  T,  H,  and  h  with  their  equivalents  in  terms  of  T,  H,  and  IT,  Equations 
(65)  to  (67),  Equations  (53)  and  (64)  become: 


V2  T  +  2-7?  +  H  =  0  (69) 
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for  all  x  C  b 


(70) 


VnT  =  h  AT 


Of  course,  Equations  (69)  and  (70)  are  the  information  equivalents  of  the 
steady  state  heat  equations  and  boundary  conditions;  however,  in  manifestly  scale 
invariant  form,  they  imply  the  existence  of  a  "temperature"  function 


T  =  T  (x,  H,  h,  AT)  (71) 

which  is  a  function  of  "position,”  and  depends  parametrically  upon  the  "heat 
sources,"  the  "heat  transfer  coefficient,"  and  system/surrounding  "temperature 
differences.”  This  "temperature"  is  in  manifestly  invariant  form  and  contains  an 
optimal  number  of  experimental  parameters  to  be  controlled  in  scale  modeling 
measurements. 

The  dimensionless  variables  associated  with  convection  and  radiation  exchange 
at  the  ship  boundaries.  Equations  (42)  through  (52),  carry  over  to  the  time- 
independent  case  without  change,  except  that  all  boundary  conditions  conform  to  a 
steady-state  by  definition. 


PHYSICAL  SCALE  MODEL  AND  DESIGN  OF  A 
MEASUREMENTS  LABORATORY 


Figures  (1)  and  (2),  respectively,  display  the  design  of  a  generic  physical  scale 
model  and  a  measurements  laboratory  suitable  for  the  study  of  IR  contrast.  The  key 
design  features  include: 


For  the  Physical  Scale  Model 


1.  Geometric  similarity  of  architecture 

2.  Geometric  similarity  of  material  thermal  properties  k(x),  p(x),  c(x) 

3.  Geometric  similarity  of  coatings  absorbtivity/emissivity  a(x)/e(x) 

4.  Electronically  controlled  thermal  sensing-source- sink  nodes  per  isothermal 
surface 


For  the  Measurements  Laboratory 

1.  Temperature  controlled  liquid  basin  and  gaseous  atmosphere 

2.  Wave  making  machine  for  scaled  waves 

3.  Rotatable  towing  bridge,  experimental  platform,  and  accessories 
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Fig.  1.  Generic  physical  scale  model 


The  design  of  the  physical  scale  model  and  measurements  laboratory  is  based 
upon  the  scaling  laws  and  dimensionless  quantities  that  transform  the  heat  equation 
and  its  boundary  conditions  to  scale  invariant  form.  In  the  most  general  time 
dependent-case,  the  scaling  laws  are  given  by  Equations  (23)  to  (26)  and  the 
dimensionless  quantities  by  Equations  (27)  to  (31). 

Consider  first  the  design  of  the  generic  physical  scale  model;  see  Figure  1. 
Geometrical  similarity  of  the  scale  model  and  the  ship  prototype  is  fundamental.  This 
is  enforced  by  requiring  distances  to  be  scaled,  so  that  there  is  a  one-to-one 
correspondence  between  geometric  points  of  model  and  ship.  The  model  is  designed 
so  that  the  dimensionless  distance,  Equation  (27),  is  the  same  at  analogous  points  of 
model  and  ship. 

In  summary,  in  laboratory  size  the  scale  model  is  to  be  a  geometrical  replicate, 
inside  and  out,  of  a  ship  prototype.  Thus,  it  is  not  adequate  to  simply  model  the 
external  surfaces  of  the  ship.  Instead,  the  interior  compartmented  architecture  is 
replicated  as  well.  Thermal  similarity  requires  that  the  thermodynamic  properties  of 
the  building  materials,  k,  p,  c,  and  e,  shall  have  the  same  spatial  dependence  in  model 
and  ship.  For  example,  the  conductivity’s  functional  dependence  upon  the  "position,” 
x,  shall  be  the  same  in  model  and  ship.  Similar  remarks  apply  to  the  density, 
specific  heat,  emissivity,  heat  sources,  coefficients  of  heat  transfer,  and  the 
temperature  differences  with  the  surroundings.  The  last  of  these  parameters  are  not 
solely  physical  properties  of  the  model,  but  obviously  depend  on  the  environment 
with  which  it  interacts;  the  interaction  must  have  the  same  functional  dependence  on 
"position”  in  model  and  ship. 

In  addition  to  geometrical  similarity  and  thermal  similarity,  the  model  is 
designed  in  a  way  that  allows  the  necessary  measurements  to  be  made  conveniently. 

A  quantity  of  fundamental  importance  in  IR  is  the  temperature,  and  it  should  be 
measured  in  various  independent  ways.  A  calibrated  IR  radiometer  (or  FLIR)  is 
installed  in  the  measurements  laboratory  for  this  purpose.  In  addition,  the  scale 
model  is  lined  with  numerous  experimental  nodes  that  serve  as  thermal 
sensing/source/sink  stations.  In  Figure  1,  the  model  is  divided  into  isothermal 
surfaces,  and  each  number  represents  an  S/S/S  station.  Each  node  carries  out 
multiple  functions  under  prearranged  computer  control.  A  thermocouple  sensing 
device  indicates  the  temperature.  A  heating  coil  serves  as  a  heat  source.  A  Peltier 
effect  junction  serves  as  a  heat  sink.  If  the  geometrical  area  of  the  experimental 
nodes  is  sufficiently  small  compared  to  the  total  external  area  of  the  model,  the 
perturbation  to  the  IR  signal  will  be  negligible,  as  required.  The  surface  of  the  model 
that  forms  the  boundary  to  the  surroundings  has  a  sufficient  number  of  nodes  attached 
so  that  each  approximately  isothermal  region  will  be  thermally  sampled.  This  also 
applies  to  the  internal  compartments;  thus  walls,  ceilings,  and  floors  are  outfitted  with 
S/S/S  nodes,  as  shown  in  Figure  1.  The  design  of  the  measurements  laboratory 
(Figure  2)  is  such  that  effects  due  to  a  ship’s  surroundings  can  be  replicated. 
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Fig.  2.  Basic  features  of  a  scale  model  measurements  laboratory. 
(Patterned  after  the  Maneuvering  and  Sea-Keeping  Facility,  CDNSWC) 


Thus  tiie  movements  of  fluids  (e.g.,  water  and  air)  past  the  ship  model  surfaces  can 
be  controlled  experimentally.  This  could  be  done  by  moving  the  model  through 
stationary  fluids  or  by  moving  the  fluids  past  a  stationary  model.  There  are  some 
experimental  advantages  in  either  of  the  measurements  methods.  But  here  a 
laboratory  design  choice  in  which  arrangements  are  made  to  move  the  model  through 
the  laboratory  past  the  relevant  fluids  is  suggested. 

A  rotatable  bridge  for  experimentation  is  suspended  over  the  liquid  basin.  A 
towing  car  that  runs  on  tracks  the  length  of  the  bridge  provides  an  experimentally 
chosen  velocity  to  the  scale  model  under  tow.  Wave  making  apparatus  along  two 
adjacent  and  orthogonal  laboratory  walls  supply  waves  of  chosen  scale  and  direction. 
The  velocity  direction  of  the  towed  model  is  also  under  computer  control  since  the 
orientation  of  the  bridge  is  variable.  Wind  effects  are  replicated  by  suspension  of  a 
mechanical  wind  machine  from  the  towing  car.  The  wind  velocity  magnitude  and 
direction  are  measurable.  The  fluids  used  in  the  model  basin  need  not  simply  be 
water  and  air.  By  using  fluids  with  differing  Re,  Pr,  and  Gr  numbers,  it  is  possible 
to  vary  the  heat  transfer  coefficients  hc  (Re,  Pr,  Gr)  to  values  convenient  for  a  given 
experiment.  For  reasons  of  feasibility,  an  appropriate  location  for  such  a  measurement 
must  be  carefully  planned.  The  temperature  of  the  fluids  may  be  controlled.  Solar 
heat  load  effects  are  replicated  by  means  of  a  carbon  arc  light  source  with  a  frequency 
spectrum  similar  to  that  of  the  sun.  The  "solar  source"  is  suspended  from  the  towing 
car  in  a  geometric  configuration  which  allows  for  variable  zenith  and  azimuth  angles, 
measured  from  the  position  of  the  scale  model.  The  intensity  of  the  "solar  source" 
may  be  varied  so  that  the  heat  load  at  the  model  surface  is  an  experimental  variable. 

A  computer  workstation  in  a  control  room  atop  the  bridge  is  used  to  "initialize"  each 
experiment  through  automated  readout  and  control  of  the  S/S/S  nodes,  to  collect  the 
data  of  each  run,  and  to  analyze  the  data  collected. 


SUMMARY,  CONCLUSIONS,  AND  SUGGESTIONS 
FOR  FURTHER  RESEARCH 

The  feasibility  of  studying  the  IR  signal  of  surface  ships  by  means  of 
measurements  upon  physical  scale  models  has  been  shown.  This  follows  from  analysis 
of  the  heat  equation  and  its  boundary  conditions  which  leads  to  the  laws  of 
thermodynamic  similitude,  i.e.,  the  scaling  laws  for  the  IR.  The  scaling  laws  have 
been  derived  with  wide  generality,  applying  to  ships  with  position  variable  thermal 
properties,  realistic  compartmented  architecture,  general  heat  loads,  convection  and 
radiation  heat  exchange  mechanisms  with  the  surroundings,  in  transient  or  steady  state 
conditions.  The  scaling  laws  are  suggestive  of  a  set  of  dimensionless  variables  which 
render  the  heat  equation  into  scale  invariant  form  whose  descriptions  therefore  apply 
to  geometrically  similar  ships  and  scale  models  of  arbitrary  size.  True  scale  modeling 
is,  therefore,  seen  as  an  activity  based  upon  experimental  arrangements  so  ordered 
that  precisely  the  same  dimensionless  magnitudes  apply  to  ship  and  model.  The 
validity  of  scale  modeling  is  ensured  when  "temperature"  measurements,  as  a  function 
of  the  dimensionless  variables,  fall  on  the  same  experimental  curve  for  both  ship  and 
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scale  model.  The  necessity  of  equating  dimensionless  variables  in  the  ship  prototype 
and  its  scale  model  replicate  leads  directly  to  the  design  features  of  both  scale  moael 
and  measurements  laboratory.  For  the  scale  model,  design  features  include  geometric 
similarity  and  thermodynamic  similarity  in  accordance  with  the  derived  scaling  laws. 
Each  important  isothermal  surface  throughout  the  model  is  attached  to  an  electron¬ 
ically  automated  sensing-source-sink  thermal  node.  Each  node  serves  the  multiple 
purpose  of  detecting  the  temperature  and  in  initializing  procedures,  either  absorbing 
or  supplying  heat,  as  requited.  For  the  measurements  laboratory,  design  features 
include  a  controlled-temperature  liquid  basin  and  gaseous  atmosphere,  scaled  wave 
making  capacity,  solar  and  wind  simulator,  a  FLIR,  an  experimental  platform,  and  a 
model  towing  capacity.  The  platform  would  include  a  computer  workstation  to 
control  the  automatic  collection  of  measured  data  and  to  carry  out  subsequent  data 
analysis. 
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